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ASYMMETRIC SIMPLE EXCLUSION PROCESS WITH 
OPEN BOUNDARIES AND KOORNWINDER 
POLYNOMIALS 


LUIGI CANTINI 

Abstract. In this paper we analyze the steady state of the Asym¬ 
metric Simple Exclusion process with open boundaries and second 
class particles by deforming it through the introduction of spec¬ 
tral parameters. The (unnormalized) probabilities of the particle 
conhgurations get promoted to Laurent polynomials in the spec¬ 
tral parameters and are constructed in terms of non-symmetric 
Koornwinder polynomials. In particular we show that the parti¬ 
tion function coincides with a symmetric Macdonald-Koornwinder 
polynomial. As an outcome we compute the steady current and 
the average density of first class particles. 
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1. Introduction 

The Asymmetric Simple Exclusion Process (ASEP) is a paradigmatic 
example of an out of equilibrium systems [36l [161 El 113 • simplest 

form the ASEP consists of particles located on the sites of a directed one 
dimensional lattice under the condition that each site can be occupied 
by at most one particle, so that the local conhgurations on each site 
can be denoted by a • for an occupied site and by a o for an empty 
site. The particles are subject to a stochastic evolution which consists 
of jumps of a particle on an empty neighboring sites with rates or 
t~2 depending on the direction of the jump. 
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Despite its apparent simplicity, this model can be successfully employed 
to describe or at least to capture the main features of very different 
physical systems. Indeed the ASEP has appeared for the hrst time 
in biology [2S] but since then it has found applications in the study 
of a wide range of physical phenomena: traffic flow, surface growth, 
sequence alignment, etc. (see |H] for a recent review of many of these 
applications). 

The ASEP can be in particular used to model the exchange of parti¬ 
cles between two reservoirs at different chemical potential im. In this 
case one considers a lattice of finite length N and the particles can be 
exchanged with a left and a right reservoir [33] • So, besides a rate for 
right and left jumps of a particle, we have also rates a, 7, respectively 
for a particle injected in or removed from the leftmost site and 5, (3 for 
a particle injected in or removed from the rightmost site. 
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Apart for its physical interest, this model has recently drawn the at¬ 
tention also of combinatorialists. Indeed the (unnormalized) stationary 
probability of the particle configurations are given by enumerations of 
certain classes of tableaux [211 El E]. 

The model we are going to consider here has as an extra feature, the 
presence of second class particles, which are confined in the finite lat¬ 
tice, i.e. that are not exchanged with the reservoirs at the boundaries. 
These second class particles can be thought as some sort of mobile im¬ 
purities and will be indicated by *. Their jump rates on empty sites. 
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or exchange rates with usual {first class) particles are given by 



1 
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We shall call this model for short “open 2ASEP”. 

The Markov chain of the 2ASEP turns out to be integrable for a 
generic choice of the parameters introduced above. Actually it is not 
difficult to show that the ASEP (in absence of second class particles) is 
algebraically equivalent to a spin | chain of XXZ type with non diago¬ 
nal boundaries. A great deal of effort has been devoted to exploit the 
integrable structure of the ASEP, either through Bethe ansatz tech¬ 
niques or separation of variables or other powerfull methods of the 
theory of integrable systems, in order to study the full spectrum of the 
ASEP Markov matrix. Reviewing all these studies goes beyond the 
scopes of the present paper, the interested reader can have a look at 
the nice introduction of [2Zj and references therein. 

A different approach, which goes under the name of Matrix Product 
Ansatz [18], has proven to be more successful in computing the prop¬ 
erties of the asymptotic steady state of the system [3|. This approach 
consists in writing the stationary probability of any particles conhgu- 
ration as a trace of products of some auxiliary (inhnite dimensional) 
matrices that satisfy a certain algebra. For concrete computations one 
has to find manageable representations of the matrix algebra. 

In absence of second class particles, Uchiyama, Sasamoto and Wa- 
daty [10] have provided a representation of the Matrix algebra in terms 
of the Askey-Wilson polynomials, which allows for the exact compu¬ 
tation at finite size of the steady state current and of its the density. 
These quantities are given by rather simple expression involving the 
moments of the Askey-Wilson measure. Combining these results of [iQ] 
with a remarkable combinatorial expression for the stationary measure 
of particle configurations in terms of staircase tableaux [T2], Corteel 
and Williams have provided a very interesting combinatorial formula 
for the moments of the Askey-Wilson polynomials [12] (see also [ID]). 

In [3D], Uchiyama has tackled the case in presence of second class 
particles, using a representation of the matrix algebra in terms of a 
t-boson. He was able to compute the so called “partition function” as 
a certain contour integral and from it to compute the phase diagram 
of the model in the thermodynamic limit. Nonetheless his finite size 
formulas are quite complicated if compared with the one obtained in 
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absence of second class particles in [10] using the Askey-Wilson poly¬ 
nomials. 

In the present paper we put forward an approach to the computation 
of the steady state of the open 2ASEP that bypasses the Matrix Prod¬ 
uct Ansatz. Such an approach has been pioneered by Di Francesco and 
Zinn-Justin in the context of the stochastic dense 0(1) loop model [19] 
and is based just on the integrability of the model. In this introduc¬ 
tion we are going to present just the rough idea behind this approach, 
whose details are discussed in the body of the paper. 

Given our choice of jump rates we know that the open 2ASEP is 
integrable, which means in particular that its Markov matrix M can 
be obtained starting from the so called scattering matrices Si. These 
matrices will be dehned in Section |2l here we only remark that they 
depend on a set of spectral parameters z = ..., zn}, where N is 

length of the lattice and the the parameter Zi is “attached” the Ath 
site. The vector gathering the stationary probabilities of each 

2ASEP conhgurations on a system of length N in the sector with m 
second class particles, can be obtained by the eigenvector of 

the scattering matrices 


5i(z)Tjv,^(z) = 
upon specialization Zi = 1 

'^N,m = 

The advantage of dealing with Tjv,m(z) is that, this vector can be nor¬ 
malized in such a way that it is a Laurent polynomial in the variables z 
and it satishes a set of exchange-reflection relations that correspond to 
a degeneration of the quantum Knizhnik-Zamolodchikov (qKZ) equa¬ 
tions of type Cn (see [2ll |3l| for some recent results on polynomial so¬ 
lutions of qKZ equations of type Cat). We shall construct the solution 
of the qKZ equations using non-symmetric Koornwinder polynomials 
and study some of their properties (notice the recent work [38] on the 
relations between Koornwinder polynomials and the XXZ spin chain). 
We shall prove that the normalization of T(sometimes called the 
“partition function”) is given by a symmetric Macdonald-Koornwinder 
polynomial for a single column partition of length N — m. From the 
weighted partition function we shall extract current and average density 
of hrst class particles. Finally, using a known integral representation 
for the Macdonald-Koornwinder polynomials we shall compute these 
quantities in the thermodynamic limit, in the case of a hxed density of 
second class particles, recovering the phase diagram of [39] . 
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The paper is organized as follows. In Section [2] we shall discnss 
the integrability of the open 2ASEP. In particular we explain how to 
make to model inhomogeneous by introducing the spectral parameters 
z and show in Theorem [1] that the inhomogeneous generalization of the 
stationary probability T(z) satishes a set of exchange-reflection rela¬ 
tions. The algebraic structures behind these equations will be discussed 
in Appendix where we also recall some basic facts about symmet¬ 
ric Macdonald-Koornwinder polynomials, non-symmetric Koornwinder 
polynomials and Askey-Wilson polynomials. In Section [3] we analyze 
the solution of the qKZ equations. We show in Theorems [21IH that the 
normalization of T(z) and the generating function for the total density 
of hrst class particles are given by symmetric Macdonald-Koornwinder 
polynomials. In Section 0] we shall discuss some recursion relations 
between solutions of the exchange-reflection equations for systems of 
different size and different values of the parameters a, b, c, d. In Section 
|5] we compute the average current and density of hrst class particles, 
we shall then take the thermodynamic limit and compute the phase 
diagram for the model. 


2. Integrability and exchange-reflection equations 

The Markov matrix which generates the time evolution of the 2ASEP 
with second class particles can be written as a sum of local terms. Let 
us associate to each site a vector space 1/ ~ with basis vectors 
{n_i,noWi} which correspond to the following local particle conhgura- 
tions 


V-i <—> o empty site 
Vq i—> * second class particle 
Vi i—> • hrst class particle 

The vector space with a basis labeled by global particle conhgurations 
has a tensor product structure = Vi ® V 2 ® ® Vjsi, with Vi ~ 

K ~ C^. The local contributions to the Markov chain are given by the 
following operators e G End(K 0 K), f(p, a) G End(K) 

1 

e = ^ (g) 0 ^|signO-i) 

ij=-l 

f (p, a) =p + a 


(1) 
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where the operators ^ End(t/) are dehned by 

The Markov matrix Ai reads 

N-l 

(2) 7W = 5^e, + fi(7,a) + fA,(/3,5). 

i=l 

where 

ej = li (g) • • • 0 lj_i 0 e 0 li +2 0 • • • 0 Iat 
fi(p, a) = li 0 • • • 0 li-i 0 f(p, a) 0 Ij+i 0 • • • 0 Iat 
Since the dynamics preserves the number of second class particles, we 
have a splitting of ifAf as a direct sum 

N 

HN = ^HN,m 

m=0 

of subspaces (sectors) Hj^^rn with hxed number m of second class par¬ 
ticles, which are preserved by the action of Ai. 

A hrst glimpse of the algebraic structure behind the Markov ma¬ 
trix M. comes from the observation that the operators ej,fi( 7 ,a) and 
fAr(/3,5) generate a representation of "Hat, the affine Hecke algebra of 

type Cat. Indeed, the operators Tq, Ti,..., Tat dehned as 

_ 1 _ 1 1 _ 1 
To = a 2 ^ ^1 

Tn = 

and for 1 < z < — 1 

Ti = ~\~ t 2 1 

satisfy the commutation relations of the generators of "Hat (see Section 

EH) 

(3) T,-T-^ = ti-t-^ 

(4) = if \^-j\>l 

(5) TiTi+iTi = Ti+iTiTi+i if z 7 ^ 0, - 1 

(6) ToTiToTi = TiToTiTo 

(7) TnTn_iTp^Tp^_i = TAr-iTArT/v-iTAr 

with ~ and tj = f for 1 < z < A^ — 1 . 

Operators written as sums of generators of Hn, like A4 are known to 
be quantum integrable [35l [20]. In what follows we shall present only 
a few notions of the theory of quantum integrable systems needed in 
order to study the properties of the stationary measure. 
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The first ingredient are the R matrices based on the baxterization of 
the Hecke algebra of type which read 

(8) Ri{z) = 1H—1- 

t'^Z — t 2 

They satisfy the braided Yang-Baxter equation (YBE) 

Ri{yz~^)Ri+i{xz~^)Ri{xy~^) = Ri+i{xy~^)Ri{xz~^)Ri+i{yz~^) 
and the so called unitarity condition 

Ri{z)Ri{z~^) = 1 . 

The second ingredient are the boundary scattering matrices. A classi- 
hcation of all the integrable boundary scattering matrices compatible 
with the matrices Ri{z) has been obtained recently |T3], here we shall 
consider the case that corresponds to the baxterization of the boundary 
matrices fi((5, a) and fiv(/ 5 , 7 ). For f 7 ^ 1 it turns out to be convenient 
to parametrize the boundary rates as 

— t~^)ab 

“ ” (a-l)(6-l)' ^ 

W/ / 1 _lx , ._i .1 

_ {t 2 -t 2 )cd t ^ - h 

^ (c-l)(d-l)’ (c-l)(d-l) 

This means that for f > 1 we can choose a, c < 0 and b,d > 1, while for 
0 < f < 1 we have a, c < 0 and 0 < b,d < 1. The boundary scattering 
matrices read 

( 10 ) A:i(^|a, b) = 1+ 7 ~^fi( 7 , a), 

[z — a) [z — b) 

(11) A'„(2|c.o!) = 1 + --AYA— 

[cz — l)(a^ — 1 ) 

They satisfy the boundary Yang-Baxter equations (BYBE), also called 
reflection equations [351 [ 8 ] 

( 12 ) R^{xy-^)K^{y)Rr{x-\-^y)K^{x) = 

Ki {x)Ri {x~^y~^)Ki {y)Ri {xy ~^), 

(13) RN-i{xy~^)Kj^{x)RN-i{xy)Kj^{y) = 

Ki^{y)RM-i{xy)Kj^{x)RN-i{xy~^), 

and the unitarity conditions 

Ki{x)Kr{x-^) = 1 , Kn{x)Kj,{x-^) = 1 . 
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Here is the explicit form of the K matrices acting respectively on the 
hrst and last site 


/ 1 0 0 

(14) Ki{z\a,b) =010 
\ 0 0 1 


z ^-1 

{z — a) {z — h) 


—1 0 —ah \ 

0 0 0 , 
1 0 a6 y 


(15) 

Kr{z\c, d) 


1 0 0 
0 1 0 
0 0 1 


l-z"^ 

(cz — l)(dz — 1) 


cd 0 1 \ 

0 0 0 
—cd 0 —1 y 


The most common way to exploit the algebraic properties of the R 
and K matrices is to use such matrices as building blocks of the so called 
double row transfer matrix, which depends on a spectral parameter y 
and which commute with the matrix M. for any values of y. The 
diagonalization of M. is then turned into the diagonalization of the 
transfer matrix [35] . 

Here instead we take a slightly different approach, which will work 
effectively for the the determination of the stationary state of A4, and 
is analogous to the one already employed in [5], for the stochastic dense 
0(1) loop model with open boundaries. Instead of using double row 
the transfer matrix, we use the so called Scattering Matrices which are 
dehned by 


(16) 5i(z) = Ri_i{ziZ~}^) ■ ■ ■ R2iziZ2^)Ri{ziZY^) 

Ki{z~^)Rl{ZiZi) ■ ■ ■ Ri-2{ZiZi-2)Ri-l{ZiZi-i) 

Ri{ZiZi+i) ■ ■ ■ RN-2{ZiZN-l)RN-l{ZiZN) 

RN{,Zi)Rf\l—\(yZi! ^jv) ■ ■ ■ Ri+li^ZijZij^2)Rii^ZiI ^j-|-i), 

where z = {zi,Z 2 , ■ ■ ■ ^zjsi}. Even though the matrices iSj(z) do not 
commute with A4, they commute among themselves 

[‘5*(z),5j(z)] = 0, 

and in the limit z —)■ 1 their eigenvectors coincide with the eigenvec¬ 
tors of M.. Indeed, let 'h(z) be a common eigenvector of iSi(z), with 
eigenvalues Aj(z) 

(17) 5i(z)T(z) = Ai(z)T(z). 

Call iSi(z) := Si{zj^i = l,Zj = z),A(z) = A{zjjLi = l,Zi = z) and 
^i(z) := = l,Zi = z), so that the specialization zj^i = 1 and 

Zi = z of eqs. dTTD reads 

5i(2:)Ti(z) = Ai(z)'f/i(z). 




L. CANTINI 


Differentiating this equation with respect to ^ and then setting z = 1 
one gets 

S'(l)>I>i(l) + S.(1)<I<'(1) = A,(1)<I<'(1) + A'(1)>I>.(1) 

but iSj(l) = 1, Aj(l) = 1 and it s simple to check that 

(18) S'(l) = . . A/1, 

t2 —t 2 

In this way we have turned the problem of the diagonalization of A4 
into the simultaneous diagonalization of the Scattering Matrices Si{z). 

The common eigenvectors of iSj(z) have nice covariance properties 
under the action of the R and K matrices. Let 'h(z) be a common 
eigenvector of Si{z) with eigenvalues Aj(z) and dehne 

Ti(z) ;= 

Ti(z) := A'i(^i)T(z), 

TAr(z) := Kn{zn)'^{'z.). 

The vectors Tj(z),\ki(z) and TAr(z) are common eigenvectors respec¬ 
tively of 

Si5j(z), So5i(z) and SAr5i(z), 
where Sj for z 7^ 0, acts on functions of z by exchanging Zi -H- 
while So, stv exchange respectively zi O zn -H- z]^^. 


In particular 



(Si5j(z)) Ti(z) = Aj(z)Ti(z) for j 

7^ z,z -M 

(19) 

(Si5j(z)) Ti(z) = Ai+i(z)Ti(z) 
(Si5i+i(z)) Ti(z) = Ai(z)Ti(z) 


(20) 

(so*Sj(z)) Ti(z) = Aj(z)Ti(z) for 

(so*Si(z)) Ti(z) = Ah^(z)Ti(z) 

3 + 1 

(21) 

isNSjiz)) TAr(z) = Aj(z)TAr(z) for 

(Sai57v(z)) T7v(z) = A)^^(z)Tai(z) 



Eqs. (119112T]) are consequences of the following commutation relations 
that are immediate to verify using the YBE, BYBE and unitarity equa¬ 
tions 

Ri{ziZ~^^)Sj{z) = (Sj5j(z)) Ri{ziz:[^^) for j 7^ z, z + 1 

( 22 ) Ri{ziz:r^^)Si{z) = {siSi+i{z)) Ri{ziZl^{) 

Ri{ziZl_^^)Si+i{z) = {siSi{z)) Ri{zizz_^^) 
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(23) 

(24) 


Ki{zi)Sj{7) = (so5j(z)) for j 
Ki{zi)Si{z) = (so*Sf^(z)) 

Kjv(2;jv)5j(z) = (sjv5j(z)) for j N 
Kn{zn)Sn{,-Z‘) = (so*S^^(z)) 

Since the R and K matrices are stochastic and preserve the subspace 
the stationary state in this sector is lifted to 4'Ar_m(z), the unique 
solution of the following equation on -ffv.m 

(25) 5i(z)4'jv,^(z) = 4'jv,^(z). 


Though such equations are still to difficult to deal with, they tell us that 
tl/ 7 v,m(z) can be normalized to be polynomial in the spectral parameters 
z. To really get a control on the solution of eqs. fl25l) we need the 
following 


Theorem 1. The solution TAr_m(z) of eqs. can he normalized in 
such a way that it satisfies the following set of exchange-reflection equa- 


tions 



(26) 

Pi ('2'i'2j_|_i ) 'h A7,m(z) 


(27) 

Kfizi)^ N,miz) 

= Sol&V,m(z) 

(28) 

KN{zN)'^N,m{'^) 

= SAr\kAr^m(z 


Proof. As we argued above, TAr,m(z) can be normalized to be a polyno¬ 
mial in the spectral parameters z with no overall factor. Call such mini¬ 
mal polynomial solution TAr_m(z). Using eqs. fll91l2ip and the uniqueness 
of the solution of eqs. fl25|) we conclude that TAr,m(z) satisfies eqs. (1261- 
[28|) . up to a multiplicative factor 

.Rj(^iA+l)^Y,m(z) fii^z'j Sj4/jv,yj(z) 

(29) Ki{zi)fit N,miz) = /o(z) SoTjv,m(z) 

A'v(^v)^N,m(z) = /Ar(z) SArTAr_m(z). 


Using the analytic structure of the right hand side of eqs. ([29]) and the 
unitarity conditions, we conclude that the functions ffiz) can take the 
form 


/i(z) 


/o(z) 

/y(z) 


(_l)«i 


( 1 _ 1 
t^Zj+l -t 2Zi 

tlzi - t~^Zi+i 




for 


(_l)KO^rfl 

(-i)""4" 


ni - azfil - bz) Y° 

V {z-a){z-b) ) 
f (z-cXz-d) Y 
“ cz)(l — dz)) 


i ^ 0 ,A^ 
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for some Ki,ei G { 0 , 1 } and where di is the degree of \I/Ar,m(z) in Zi. 
Setting Zi = Zi+i in the hrst of the eqs. fl291) leads to 

^N,m{'^)zi=Zi+i — (~ 1 ) N,m{'^)zi=Zi+i 

which forces the k* = 0 , otherwise one would hnd that ^Ar^m(z) is 
divisible by Zi — Zi+i. In the same way setting 2:1 = ±1 in the second 
two or zn = ±1 in third of eqs. fl 2 ^ one get kq = kn = 0 and moreover 
that di and d^ must be even integers. 

The last step of the proof consist in combining eqs. fl2^ and eqs. fl2^ . 
This leads to e, = 0 and di = d^ = d, therefore it is immediate to 
conclude that the Laurent polynomial 

N 

'hjv,m(z) := ^N,m{^)zi=Zi+, 

i=l 

satisfies eqs. fl261l28|) . □ 

A simple consequence of Theorem [1] is the following uniqueness result 

Theorem 2. The system of equations admits a Laurent poly¬ 

nomial solution, unique up to a multiplication of a function invariant 
under the action o/Sj. 

3. Analysis of the solution of the exchange-reflection 

EQUATIONS 

In order to analyze the solution of the exchange-reflection equations 
fl261l28p we need to introduce some algebraic background. 

3.1. AfRne Hecke algebra of type Cn and Noumi representa¬ 
tion. The affine Weyl group Wat of type Cat is the group generated by 
elements Sq, Si,..., sat subject to the relations 

(30) = 1 

(31) SiSj = SjSi if |i-j|>l 

(32) SjSj+iSi = Si+iSiSi+i if i 7 ^ 0, iV - 1 

(33) SoSiSqSi = SiSqSiSq 

(34) SATSAT-lSArSAT-l = SAT-lSArSAT-lSAT 

The hnite Weyl group of type Cn is the subgroup of Wn generated 
by elements Si,... s^. The affine Weyl group Wat has a faithful action, 
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depending on a parameter g, on , 2 ;^], the space of Laurent 

Polynomials in N variables 

(35) Zi+i, ...,zn) = f{zi,Zi+i, Zi,..., zn) 

(36) So/( 2 :i, ...,Zi, Zi+i, ...,zn) = ■ ■ ■, Zi, Zi+i,..., Zn) 

(37) SNfizi, ...,Zi, Zi+i, ...,Zn) = fizi, . . . , Zi, Zi+i, ..., z^^) 

It has also actions on , parametrized by the level £ G Z, 

(38) • • • •) D 25 D 2 - 1 -I 5 • • • ; Djv} • • • f O^i, • • • , Djv} 

(39) * * * ; C^i-^i, • • • ; Djv} 0^1, • • • 5 O^i, * * * ; 

(40) SAr{ai, ttj+i, . . . , Oat} = {«!, . . . , Q!i, Oj+i, . . . , —ttAr}. 

Notice that the action of on , z^] can be read off from its 

action on by looking at monomials, namely for w G 


wz = z 


We dehne two partial orders on Z^. The hrst is the usual dominance 
order-, a < /? if for 1 < j < iV we have ~ f^i) — 0 - The second 

order, a ^ (3, is dehned as follows. Call a'^ the unique element in W^a 
such that is a partition i.e. such that cxj > Oj+i > 0. Then we say 
that a ^ (3 ii < (3^ OT in case a'^ = (3^ ii a < (3. 

The Affine Hecke algebra "Hat of type Cn is a deformation of the 
group algebra of Wn, which depends on three parameters to,tN,t and 
is generated by elements To,Ti,... ,Ti\f subject to the commutations 
relations 

(41) t,-T-^ = 4- 

(42) TiTj = TjTi if |i - jj > 1 

(43) TiTi+iTi = Ti+iTiTi+i if i 7 ^ 0, - 1 

(44) ToTiToTi = TiTqTiTo 

(45) = Tn-iT^Tn-iTn 

with ti = £2 = ■■■ = tjsf -1 = t. 

The hnite Hecke algebra of type Cn is a sub-algebra of Hn 
generated by elements Ti,... ,Tn-i,Tn. It is well known [23] that a 
basis of "Hat is parametrized by elements of Wn 

T = T- T- T- 

where w = Si^^Si^ ■ ■ ■ Si^ is a reduced expression of w G Wat- 
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An important commutative sub-algebra yN is generated by Lusztig 
elements ..., [28] 

(46) y = {T,... T^_i)(Tjv ... To)(Tf 1 ... Tr\). 


In [3T], Noumi introduced a representation of "Hjv depending on 6 
parameters a, 6, c, d, t, g, acting on ..., the space of Laurent 

polynomials in N variables z = 


(47) 

(48) 

(49) 


Ti=U - {t^Zi - t ^Zi+i) di 
-1 {zi - a){zi - b) 


To = to 


f ^ 


do 


Zl 


Tm — t 


-1 {czn - l){dZN - 1) 


N 


f ^ 


Zn 


dN, 


where = —q ^ah tjsi = —cd and the hnite difference operator di, Oq 
and are dehned by 


di = 


do = 


So 


- Zi+i ' 


-1 ’ 

Zl - qz^ 


d]\f — 


Sn 


Zn - ’ 


3.2. Analysis of eqs. (l26tl28h . Let us write the wave function 4/ N,m{'^) 
in the basis w = wi... wn, uii = o, *, • 

^'iv,m(z) = V'w(z) w, 

wGQ(W,m) 

where the sum runs on Q{N,m), the set of 2ASEP conhgurations on 
a strip of length N in the presence of m second class particles. We 
can identify a word in o, *, • with an element of through the rules 
o = —1,* = 0, • = +1, hence the set Q{N,m) inherits an action of 
and of Wat at level 0 and the partial orders <, ^. 

Once written in components, the exchange-reflection eqs. (12611251) are 
expressed in a nice compact form using the generators of the affine 
Hecke algebra "Hat in the Noumi representation at g = 1 

(50) iiu, = ti^fi'ip^ if 

(51) 'ipsiw = tifi'iljy, if 

Recall that for 1 < ^ — 1, f* = f, while Iq = —ah and In = —cd. 

In words, eq. fISOD tells us that if Wi = tCj+i then is symmetric 
under exchange Zi -H- Zj+i, while if tco = * or tCAf = *, then - 0 ^ is re¬ 
spectively symmetric under inversions Zi Zn - Eq.([5lD 


W = SiW 

j i > 0 & w < SiW 
[ i = 0 & w > sow 


^Here we shall just consider a, b, c, d,t,q € C. 
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tells US how to exchange two different neighboring particles 


o • • • 


• • • o 


thfi 

* • • • - ^ • • • * 


■ ■ • * o • • • 


• • • o * • • • 


and how to inject/remove particles from the boundaries. 




^nTn 


-)■ • • • O 


3.3. Reference component. Let us consider a system of length N = 
k + m, in the sector with m second class particles, and analyze the 
component associated to the state 


, m 


= {o • • • o * • • • *} 



It follows from the eas.flSOllSTj). that (z) is preserved by the Lusztig 
operators f/. Therefore (see AppendixEl) we can identify '0ofe*^(z) with 
the non-symmetric Koornwinder 

(52) V’ofc*-(z) = E^(k,m)iz) 

associated to the string 


/i(fc, m) = { 


This fixes completely the normalization of T 7 v,m(z). 

The non-symmetric Koornwinder polynomials depend on the further 
parameter q which appears in the general Noumi representation. Hence 
a priory we should specify that E^(k^rn){z) is taken for g = 1, but 
actually we will see in a moment that E^(^k,m){z) doesn’t depend on q. 

Moreover, while the monomials expansion of any non-symmetric 
Koornwinder polynomial reads as in eq. fllOip . the expansion of E^(^k,m){'^) 
is much simpler, it contains only negative powers of Zi for 1 < i < k. 

Either from the exchange eq. (I50|) applied to V’o''*"* (z) or from the gen¬ 
eral theory of non-symmetric Koornwinder polynomials [32] we know 
that E^(fc^m)(z) is separately symmetric in the first k and last m spectral 
parameters, so that we can express it as a sum of elementary symmetric 
polynomials in z variables 


E, 






ek-i{z^ \ 




i=0 


(53) 
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with = 1. Using the symmetry properties of we can 

rewrite its dehning equation (I102p as 

(54) TfcTfc+i... T/v-iT/vT;v_i... TiToU^(fc^m) (z) = E^(^k,m) (z) 

Pictorially this equation has a clear meaning: the action of Tq injects 
a particle on the left, then Tat-i ... Ti moves this particle to the far 
right, Ttv removes the particle from the right and finally TkTk+i ■ ■ ■ Uv-i 
brings back the empty state in position k. 

Equation fl54p allows to determine in a recursive way. It turns 

out that they do not depend on k and have a rather simple dependence 
on m. This is the content of the following 

Proposition 1 . The non symmetric KoornwinderpolynomialE^(^k,m){'^) 

does not depend on q and reads 


(55) E^(fc,„.)(z) = I^H{w;a,b,rc,t^d)ll 

i=i 




where the integration is around zero and H{w; a, b, c, d) is a formal 
power series in w 


(56) 


with 


H{w; a, b, c, d) = hn{a, b, c, d)w^ 

3=0 


(57) hn{a,b,c,d) = 
1 


{abed; f). 


E 


0<ij,k<n ^ 


n 


t(2)+(0a*h^'(-c)"-^(-d)*+^+^ 


J t 


where we have used the usual q-multinomial with base t 


n 

.hi, 


{t-t), 


n—1 


= ^(1 - 


Proof. Let us rewrite eq. fl53l) as a contour integral around the origin 


E, 




Z = 




2=1 


(58) 
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with Hk,m{w) = Yl'j=o ■ Using the following relations 

{w + a) {w + h) 


To{w ^ + zi)=q Uq M (w + Zi) 


w 


Tj{w + Zj){w + Zj^^) =t ^{w + Zj ){w + Zj+i) 

Tj{w-^ + 2 : 7 + 1 ) = t-^tvj-^ + zj^) 

Tj{w + Zj) = {tw + Zj^i) 

^ , , -if , (1 + wc)(l + wd) 

Tn{w + Zn) =tj^^ [{w ^ + Zjf^) - 


w 


in eq. 


we get 
h m f dw 

tlt^ ■ 


Hk 


w 


(a + w)(b + w) 


k-l 


w 


dw , ,, {1 + t'^cw){l + t"^dw) 


<'N 


^ + Zj^] 

i=i 

fc-i 


27ri 




W 


If:*"" + T‘) 

i=i 


The fact that this equation must be true for any choice of ..., Zk-i 
means that we can substitute with any polynomial 

p{w~^) of degree fc — 1, getting 


(69) 


(p —Hk,m{w) \p{w 7 - 

io 271? V w 


From the previous equation we can already make a couple of conclu¬ 
sions. 

1) Since in eq. (l5^ the parameter q (which was present through to = 
—q~^ab) has disappeared we have that E^(^k,m){'^) does not depend on 

Q- 

2) The dependence on k is only in the degree of the arbitrary poly¬ 
nomial p{w), which means that the coefficients actually do not 

depend on /c, therefore we can suppress the label k and we think of 
Hk,m{w) = Hm{w) as a formal inhnite series. 

3) The parameter m appears just by multiplying both c and d by T" 
and we can write the solution of eq. flS^ as 

Hmiw] a, 6, c,d) = H(w; a, 6, t^c, d^d) 
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It remains only to use eq. (j5^ to determine the coefficients /ij(a, b, c, d). 
By setting p{w) = into eq. flS^ we get following three-terms re¬ 
cursion relation 

(60) — c~^d~^)hn+ 

(t^ ^(a + b) — {c ^ + d ^))hn-i + ^ — l)hn -2 = 0. 

with initial conditions ho = 1. It is then not difficult to check that the 
expression in eq. (l57|) satishes this recursion relation. □ 

An immediate consequence of Proposition [T] is that 

(61) ij^Qk ^m(z^ b^ Cj d) — ^m—i (zjY; b^ tc^ td^. 

where zy means that the variable Zj is absent. 

The recursion relation for the coefficients hj eq. 0601) has some invari¬ 
ance properties. It is invariant under a O c~^,b -H- d~^,t ■(-)■ from 
which (writing explicitly the dependence of h„ on t as h„(a, 6 , c, (i|f)) 
we get the duality relation 

(62) hn{a, b, c, d\t) = hn(c~^, d~^, a~^, b~^\t~^). 

It is also invariant under 

CL —^ X (X, b —X b, c —^ Ac, d —A d, hn —^ A 

which means that h„(a, b, c, d) is homogeneous in a~^, b~^, c, d of degree 
n and we have 

(63) H{X~^w; X~^a, X~^b, Ac, Ad) = H{w; a, b, c, d). 

Using repeatedly the relation -)- zj^) = 

we obtain easily the components of all the conhgurations without hrst 
class particles. For J = {ji, • • •, jfc}, 1 < ji < 32 < ■ ■ ■ 3 k < N call 
o(J)jv the conhguration with empty space in positions J and second 
class particles in all the other sites. Then we have 

k 

For conhgurations of the kind *"* 0 ^, combining the previous equation 
with eq. fl63|l we obtain the analog of eq. fl 6 T|) 

'ip^m^k (z, u, b^ c, d) ^ Lp^m—iQk (zj,f CL^t b^ d'^. 


(64) 
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Analogous integral formulas can be written for the components of the 
conhgurations without empty sites. For later reference we write explic¬ 
itly hi 


(65) 


hi (a, b, c, d) 


a + b — c ^ — d ^ 
c~^d~^ — ab 


Another interesting remark is that from a close inspection it turns out 
that the three-terms recursion fIHUD is nothing else than the recursion 
relation for Al-Salam Chihara polynomials in disguise (see eq. fllldp i 
therefore we have also 


( 66 ) 

hn{a,b, c, d) 


(-l)Vt^f A 
[abed] t)n y2 



3.4. Other components and normalization. Once the component 
is known, all the others are obtained from it by action of the 
generators % using the exchange equations fl5T|) . Well, actually it is not 
necessary to use the full algebra being enough. Moreover the 

state is preserved by the generators T* for 1 < i ^ N — m < N. 

In order to see how to remove this redundancy and to provide a compact 
formula for any component we need to concentrate for a moment 

on the Weyl group and its action on Q{N,m). The state w{N — 
m, m) is preserved by J := Sj^^rn x hVmi which is a parabolic subgroup 
of i.e. a subgroup generated by a subset of the generators of W^. 
In the present case I is generated by Sj for — m<iV. 

The states in Q{N,m) are in bijection with elements in W^/I. On 
the other hand, from the basic properties of parabolic subgroups of a 
Coxeter group (see [23] Chapter 1) we know that the set dehned 

>V^ = {g e >V^|£(gs) > £(g),Vs e /}. 

is in bijection with classes in the quotient W^/I, or equivalently, that 
any classes in [g] G W^^// contains a unique shortest representative 
h G and therefore any g G can be written in a unique way 
as a product g = hs with h G and s G /. Moreover, any reduced 
expression of g is the product of a reduced expression of h times a 
reduced expression of s, which in particular implies that at the level of 
the Hecke algebra, we have a decomposition 

(67) Tg = ThT, 
with h G and s G /, and 

( 68 ) xm = X(TMT.). 
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The bijection between and Q{N,m) is now obvious 

h ^ ho^ 

Coming back to the components we easily check that for h G we 
have 

(69) = x(Th)Th'0o'=*-(z)- 

It is now easy to prove the following theorem that identify the nor¬ 
malization of \h 7 v^m(z) (sometimes called also the partition function) 
with a Macdonald-Koornwinder polynomial [26] 

Theorem 3. Let 

(70) ^jv,m(z) := ^ V’u,(z) 

w£Q(N,m) 

Then we have 


(71) ■^_/V,m(z) Pri{N,m) ■ 

where Pri(N,m){'^) is the Macdonald-Koornwinder polynomial [26| asso¬ 
ciated to the partition 


t]{N, m) 



N—m m 


Proof. From [32] (see Theorem [6] in Appendix |A]) we know that 

gGWjO 

Now rewrite the sum over g G as a double sum over h G and 
over s G /, using eqs. fl67ll68p 

-P)7(W,m)(z) OC EE x(Th)x(T,)ThT,i/.„.,,„(z). 
hew^ sG/ 

For s G / one has that Tgi/’ qA; (z) OC -0 qA; (z) and hence 

sei 

Therefore we obtain 

Prj{N,m){z) OC ^ x(Th)Thl/>o'=*- (z) = ^ V’«)(z). 

hGW^ w£Q{N,m) 

In order to £x the proportionality constant we recall that Pri(N,m){'^) is 
normalized in such a way that the monomial coefficient 1. 

On the other side, in the sum dehning such monomial comes 
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only from the component and has coefficient 1. Therefore we 

conclnde that the proportionality factor is 1. □ 

Actnally we can do even better and consider the generating fnnction 
for the nnmber of particles in the system. Calling •{w) the nnmber of 
particles present in the state w, what we wonld like to consider is the 
weighted partition fnnction 

(72) ZJv(^;z):= 

w£Q{N,m) 

In order to state onr resnlt it is convenient to switch to a notation 
in which the dependence on the parameters a, b, c, d is explicit writing 
ZN,m{'2‘] a, b, c, d) for the normalization. Then we have 

Theorem 4. The generating function for the number of particles in 
the system is given by 

(73) z; a, b, c, d) = a^, 6^, c^, d^). 

with a^ = f a,h^ = f b,c^ = ^~^c, d^ = ^~^d. 

This theorem is jnst an immediate corollary of the following 


Proposition 2. For any state w G Q{N,m) we have 
(74) c, d) = a^, 6^, c^, d^). 


Proof. For w = eg. (17^ is eqnivalent to eg. (1631) . Then we claim 
that if eg. (1601) is trne for some w then it is trne for all the states with 
the same nnmber of hrst particles. Indeed if •{w) = •{w'), then 
is obtained from throngh the action of the sole operators Tj with 
1 < i < iV — 1, which are homogeneous in the variables z. 

The last step is an induction on the number of hrst class particles. 
Assume eg. (l60|) is true for all w such that •{w) < k. Then consider a 
state w with •{w) = n, of the form w = wo. Writing as 

V'^o(z) = 

the relation fl60|) translates into 

_ 1 

Acting with on the component 'ipwoi'^) we get 

ipw.{z) = {zn - C~^ - d~^)(j)Z^\z) - c~^d~^(j)^^\z). 

Using this expression and egs. dT^ we easily verify that ZwmiZj satishes 
eg. (l60|) . □ 
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A more explicit formula for the weighted partition function ZN,m{e-,z), 
which is well suited for asymptotic analysis, is obtained in terms of an 
integral representation of the Macdonald-Koornwinder polynomials [30] 
(see Appendix IA.2|) . Assuming t < 1 we can write 


(76) ZN,mi^ ■,z)=r^ {a^,b^,c^,d^\t)x 




N—m 


dx 


Anix 


-n(z , x)w{x] ag, cg, d^\t)pm{x] ag, 6^, c^, d^\t) 


where Pm{x] a, b, c, d\t) is the m-th Askey-Wilson polynomial of base t 


m 


the variable —, w{x; a, b, c, d\t) is the Askey-Wilson kernel 


(77) w{x; a, b, c, d\t) = 


{x^,x 2;t)^ 


{ax, ax~^, bx, bx~^, cx, cx~^, dx, dx~^-, t)oo ’ 

while n(z,x) = ni<i<Ar(^* + — x — x~^) and the normalization 

constant is given by 

{abcdt^’^] t)oo 


(78) rm{a,b,c,d\t) = 


{t^+^,abt^, act™, adt^, bct^, 6dt™, cdt™; t). 


The contour of integration C encircles the poles in a(t^,b^t^ ,C(t^,d^t^ 
{k G Z+) and excludes all the others. 

In the homogeneous limit, our result for Z 7 v,m(^;z) looks quite dif¬ 
ferent and simpler than the analogous result in [3D], where Uchiyama 
computes the grand partition function. A more direct comparison 
can be done for m = 0 with the expression for ZM^rn{i]^) presented 
in [ini Eq.(6.2)] (their parameters a,b,c,d,q are related to ours by 
q = t~^,a t-A —a~^,d y-)- —d~^,b -f-)- —c“^). Since po{x\t) = 1, we see 
that the integral in eq. fl76|) coincides with the one in [lO] Eq.(6.2)], 
nonetheless the two formulas differ by a prefactor which depends on ^ 
(and therefore cannot be accounted by a different normalization of the 
wave function), but not on the size N of the system. This gives rise at 
hnite size to different predictions for the average occupation number of 
hrst class particles. 


4. Recursion relations 

In this Section we will discuss some relations between solutions of the 
exchange-reflection equations for systems of different size and different 
values of the parameters a, b, c, d. This will allow us to derive some 
contiguous relations for Askey-Wilson polynomials. Whenever we shall 
compare solutions corresponding to different values of the parameter 
a, b, c, d we shall keep explicit the dependence on these parameters. 








ASEP WITH OPEN BOUNDARIES AND KOORNWINDER POLYNOMIALS 21 


Theorem 5. The following three kinds of recursive relations holds 


(79) 


'il}u,*{2;a,b,c,d) = tljy,{z^]a,b,tc,td) 
a, 6, c, d) = tb, c, d). 


(80) 


V’^o(z) + cd'4!w,{z) 

ab'fowi'z) + V'.«,(z) 


(1 - czn ){1 - dzN) 
Zn 

(a - zi){b - zi) 


'*/’ui(ZjY) 


2:1 




(81) 


(82) 


V'«,o(z;a, 6, c,(i)|^^=c-i = KR{c)fjw{zfj;a,b,tc,d) 
^p^oiz;a,b,c,d)\^^=d-^ = KR{d)'ipy,{z{^-,a,b,c,td) 
V>^,(z;a,6,c,rf)|2j^=c-i = -c~^d~^KR{c)'ijj^{zf^; a,b,tc, d) 
^p^,{z;a,b,c,d)\^^=d-^ = -c~^d~^KR{d)'4)^{z}^-,a,b,c,td) 

fjowia, b, c, d; z)\^^=a = KL{a)ii^{ta, b, c, d; Z|) 
fjowia.b, c,rf;z)| 2^=6 = KL{b)'ijj^{a,tb, c, d; zj) 
'ip,^{a,b,c,d]z)\^-^=a = -abKL{a)'ilJw{ta,b,c,d]Zi) 
'f,uj{a,b,c,d-,z)\^^=b = -abKL{b)'il)^{a,tb, c, d] z{) 


where Z| means that the variable zi is absent and 


(83) 


Kr{x) 

Kl{x) 


{1 — axt^){l — bxt^)cdx ^ 
1 — abcdt"^^ 

(1 — cxt'^){l — dxt'^)x~^ 

1 — abcdt'^'^ 


Proof. First we claim that if eqs. (179118^ hold for some w G Q{N,m) 
then they hold for any conhgurations in Q{N,m). Indeed, for the 
equations in which the particle conhguration at site N is hxed, one can 
modify w by acting on both sides of the equalities either with operators 
Ti with 0 < z < — 1, while for the equations in which it is the hrst site 

conhguration to be hxed, one can modify w by acting with operators 
Ti with 1 < z < iV. 

We shall use the previous remark to prove eqs. fl79ll80p . We could use 
it also to prove eqs. fl81ll82p but we prefer to adopt a diherent strategy. 

The hrst equation of fl79p . in the case w = coincides with 

eq. dhTp . while the second equation of d79P for w = *”^ 0 ^ coincides with 
eq. dMp . 
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In order to prove the second of eqs. (1801) for w = we use the 

integral expression of b, c, d; z) given by eq. fISSD 


1 

{ab + tQTo)'ijjok^m{z) = 

{ab + f|To) ^F^w) n («;-' + zj^) 


(a — Zi){b — zi) f dw 


i=i 


Zl 




(a - zi){b - Zl) 


Zl 


(^i) 


where in the third equality we have used the identity 


/ , _iN (a — zi)(b — Zl) 

{ab + tQTo){w ^ + -^1 ^) =-. 

Zl 

In a similar way one proves the hrst of eqs. flSOp . 

For the proof of eqs. flHTj) we adopt a different strategy: we show that 
the l.h.s. satishes the same exchange-reflection equations as the r.h.s. 
and therefore by Theorem [2] they must be proportional. 

Lets look at the hrst of these equations, since we are comparing 
systems of different length we must distinguish between operators rep¬ 
resenting the Affine Hecke algebra of different rank and different pa¬ 
rameters therefore just for the present proof we shall use an heavier 
notations. For the case of rank N with parameter c highlighted we 
shall write the generators as T^^\c). In terms of divided difference 
operators we have of course T^^\c) = T^^~^\c') for i < N — 2 and any 
c, c', therefore both sides of the hrst of eqs. flHT]) satishes eqs. fl50ll5ip for 
i < N — 2 and Tj^^\tc). It remains to prove that the r.h.s. satishes 
eqs. d^TP for i = A^ — 1. If w = w* then 'ip^o does not depend on zp^^i so 
in particular it is invariant under z^^i •H- z]^^_-^^. li w = wo then ^jJu,oo 
is symmetric in zn_i,zn and we have 


(84) 'i/'*.o(z) = V’t5oo(z) 

It is straightforward to verify that for any symmetric function H{x,y) 
one has the following identity 

= (T^Z~i'\tc)y' H{zn-i,c-^). 
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This means that after specialization zn = c ^ in eq. IjMjl we get 

'0'i5»o(z) |2^=c-l (z) |2^=c-l • 


It remains to hx the proportionality constant Kr{c). This is given by 
the coefficient of the monomial 11^=1 ^7^ a, c, 

since the same monomial has coefficient 1 in (z^; a, b, tc, d). We 

get 


Kr{c) = c + rhi(a, 6, rc, t^d\t) = 


T + acf”^)(—1 + bct'^)d 


—1 + abcdf^^ 

In a similar way one proves the identities for the specializations zi = 
a, b. □ 


Since '0^o(z; a, b, c, d) and '0u„(z; a, 6, c, d) are polynomials of degree 
1 in zjq and z]^^ respectively, they are recursively determined by the 
two equations liiii. At the level of partition function the recursions 
fl79]l8B imply 


(85) ^JV,m(z, U, 5, C, d) 1 ^7V—l,m—l(Zjy,Cl,5,fC,td)T 
(l-acT”)(l-6cf”^)(l-dc) 


c(l — 


l,m —1 (zjyU, 5, fc, d). 


By using the integral formula (|7T1) in the previous equation, we derive a 
contiguous relation for Askey-Wilson polynomials which is not difficult 
to prove by direct mean^ 


(86) (1 — cdt'^)Pm{x] a, b, c, d|t) — (1 — cd)Pm{x; a, b, tc, d\t) 

— cd(l — f™')(l — abt"^~^){z + z~^ — d — d~^)Pm.-i{x-, a, b, tc, td\t) = 0. 


On the other hand eq. flHSl) and its analog for zn = d“7 determine 
ZN,mi'^',a-,b,c,d) uniquely. Therefore assuming eq. fl86l) . we could ob¬ 
tain a proof of the integral formula flTTl) . without employing Mimachi’s 
formula. 

Another consequence of the recursions fl81ll82|) concerns the coeffi¬ 
cients hn in the expansion of Indeed by applying eqs. fl81ll82p to 

the case w = we get 


hn+i{a, b, c, d) -f chn{a, b, c, d) 


(1 — ac)(l — bc)d 
(1 — abed) 


, hn{a, b, tc, d) 


hn+i{a,b,c,d) + a ^hn{a,b,c,d) = 


(1 — ac)(l — ad)a ^ 
1 — abed 


-hn(ta, b, c, d). 


^Several contiguous/difference relations for Askey-Wilson polynomials emerging 
from our discussion are presented in Appendix IA.3I 
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Actually these two equations are equivalent thanks to the duality rela¬ 
tion flB^ . The previous two equations can be written as contiguous or 
difference equations for Al-Salam Chihara polynomials 

(87) aQn+i{z; a, b\t) - (1 - abt'^)Qn{z; a, b\t) 

— a{z + z~^ — a — a~^)Qn{z-, ta, b\t) = 0, 


(88) Qn+iiz; a, b\t) - z{l - abt'^)Qn{z; a, b\t) 


Qn{t^ z;t^a,t^b\t) = 0. 


Both equations can be lifted to contiguous or difference relations for 
Askey-Wilson polynomials (see Appendix lA.Sp . 


5. Density and current 

In order to compute physical quantities, like average density of first 
class particles or current, we need to specialize the spectral parameters 
z = 1. While for the density this is straightforward 

(89) 

for the steady current {JN,m) we need to make an intermediate step 
and define an inhomogeneous version of that quantity 

(90) aV’ot«(z) - 7i/^,^(z), 

^Y,m(z) 

w&Q(N,m) 

which under specialization z = 1 reduces to the homogeneous case. 


{JN,m) '7Af,m(l)- 


Using the fact that a/y = —ab and the first of eqs. flHOjl . we get an 
explicit formula for JN,m{'^) 


(91) 


'7/V,m(z) 



Zl){b - Zi) ZN-l,m{z \ Zi) 

Zl ZN^mi,'^') 


After specialization z = 1, we obtain 






l,m(l) 

^V,m(l) 


(92) 
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5.1. Thermodynamic limit. Since the stationary current and the av¬ 
erage hrst particle density have been expressed in terms of the weighted 
partition function, in order to determine the behavior of such quanti¬ 
ties in the thermodynamic limit we need to work out the asymptotic 
behavior of 1) when N and m tend to inhnity at hxed density 

of second class particles = m/N. 

If we assume f < 1, we can use eq. fITbD . In that equation the prefactor 
cg, for m —)■ cxD goes to a constant n, therefore we just 
need to consider the integral, which in the homogeneous specialization 
reads 


(93) / da; (x^,x ^-x-x ag, cg, dg|t) 

Jc Airix (a^x, b^x~^, c^x, c^x~^, d^x, d^x~^;t)oo 

For the large N, m we can use the asymptotic formula eq. flll2p for 
Askey-Wilson polynomials and we arrive at 

(94) 

1) ke-"' 


dx (x ^;t)oo(^ + ^ ^—x — x ^)^x ^ 
c 47rix (a^x~^, b^x~^, c^x~^, d^x~^; t)oo 


For ^ 1 the function 


f(x) = log(^ + ^ ^ -x-x log(x) 

has a saddle point at 

(95) = xo + 0((^ - 1)), Xo= <-1. 

-1 + p* 

Now recall that for 0 < t < 1 the parameters a, b, c, d span the range 
a, c < 0 and 0 < 5, d < 1. When xq < a, c < 0 we can deform the 
contour of integration to the steepest descent path that passes through 
x(p*;^) and we easily hnd 


(96) 

^V,m(^^; 1) — 


k^N-m 1 

■y/87r/"(x^)iV (a5X^\6gx^\cgx^\d5X^^;f)oo 


When at least one among a and c is less then xq then the dominating 
contribution comes from the pole around if a < c or around if 
c < a, and we get 


( 97 ) 


-s— ^ ^ ^ a < c, xo 


2a^(t,b^a^ ^;t)c 


fcg^-’"(c-ht)oo(g+r^-cg-c-b^c-’ 

d'}oo 


C < a, Xq 
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Let US now call Jp^ and the steady state current and density of 
first class particle in the thermodynamic limit. Combining eqs. fl89ll92p 
with eqs. (19611971) we obtain 


(98) Jp^ 


{t2 -t 2){l - pI) 
4 


Pp 


I- p* 
2 


for xo < a, c 


(99) 

( 100 ) 


■Jp, — 


Jp* — 


, _ 1 1 . 
a[t 2 — t 2 ) 

(1-a)^ ’ 

c{t~ 

“a 


'2 — t2' 


Pa, = 


a — 1 


P* 


Pp. = 


Upon redehnition q = t a -H- —a d —d ^,60 
are in agreement with 


for a < Xq,c 

for c < xo,a 
—c~^ our results 


6. Conclusion 

In this paper we have analyzed the steady state of the Asymmet¬ 
ric Simple Exclusion process with open boundaries and second class 
particles by deforming it through the introduction of spectral param¬ 
eters, in a way which is dictated by the integrable structure of the 
model. The (unnormalized) probabilities of the particle conhgurations 
get promoted to Laurent polynomials in the spectral parameters, that 
are constructed in terms of non-symmetric Koornwinder polynomials. 
Moreover we have shown that the partition function coincides with a 
symmetric Macdonald-Koornwinder polynomial. As an outcome we 
have computed the steady current and the average density of hrst class 
particles. It is interesting to pursue further the analysis of the inhomo¬ 
geneous model and compute other quantities like correlation functions. 

In a recent preprint |T3] Corteel and Williams have uncovered a dif¬ 
ferent connection between open 2ASEP with open boundaries a second 
class particles and the theory of Koornwinder polynomials. It would be 
extremely interesting to investigate the relation of their findings with 
the algebraic structure presented in our paper. 

It is clear that our approach can be applied to other stochastic in¬ 
teracting particle systems as soon as one recognizes that the Markov 
matrix is integrable. The most straightforward generalization of the 
present paper would consist in a system with more than two species of 
particles and exchange rates that arise from a higher rank quotient of 
the affine Hecke algebra of type Cn- A more interesting class of systems 
that can be treated with our approach | 3 ], arise from the classihcation 
of the stochastic boundary scattering matrices of Crampe et ah [15] . 
One could also consider systems with periodic boundary conditions. 
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While for just a single class of particles the steady state is simply uni¬ 
form, this is no longer true if one allows for more classes of particles, 
and the nontrivial steady state can be analyzed using the approach 
presented here [3 [ 61 - 

Appendix A. Koornwinder, Macdonald-Koornwinder and 

Askey-Wilson polynomials 

In this Appendix we shall the definitions of Koornwinder, Macdonald- 
Koornwinder and Askey-Wilson polynomials, and a few of their prop¬ 
erties we needed in the paper. 

A.l. Nonsymmetric Koornwinder polynomials. In order to in¬ 
troduce the Non-symmetric Koornwinder polynomials one considers 
the action, in the Noumi representation, of the commutative sub-algebra 
yjM on ..., It turns out that such action is triangular with 

respect to the order ^ on the monomials 1371 

(101) Fiz" — w,((i)z“ + ^ 

Z“ 

in order to write the eigenvalues 0Ji{a) we need to introduce Wa, which 
is the minimal length element in such that WaC(~^ = a and p = 
{1, 2,..., iV} e Z. Then we can write 

a;.(a) = i^q-^ahcd 

The common eigenfunctions of are the non-symmetric Koornwinder 
polynomials Ea{z) [32] 

E„(z) = z" -h ^ cpz^, 

(102) z“ 

YiEa{z) = Ui{a)Ea{z). 


A. 2 . Macdonald-Koornwinder polynomials. For a partition A G 
Z^, call Hx the subspace of C[z ^,..., z^] generated by the polynomials 
Ea{z) such that = A is preserved by the action of Kn and for 
generic values of the parameters a, b, c, d, q, t it forms an irreducible 
representation. The subspace ^, invariant under the action of is 
one-dimensional and is generated by Px{z), the symmetric Macdonald- 
Koornwinder polynomial in N variables z, associated to the partition 
A [261132]. Before providing a dehnition and discussion of Macdonald- 
Koornwinder polynomials let us state one of the main results we need. 
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Theorem 6 ([32j, Corollary 6.6). Let x ■ hLw C be the one dimen¬ 
sional representation ofT-L^ defined on the generators by 

x(Ti) = tl 

Then we have 

(103) Pa+(z) oc ^ x(T^)T^E„(z) 

w&W% 

For a fixed N, the Macdonald-Koornwinder polynomials are Laurent 
polynomials in N variables, invariant under W^, which are eigenvectors 
of the following q-difference operator 

N 

P,, = - 1) + - 1) 

i=l 

where is the Fth q-shift operator 

Tq,zJ{zi, ...,Zi,...,ZN) = f{zi ,..., qzi, ...,Zn) 

and 

(1 - az){l - hz){l - c^)(l - dz) A (1 “ - ^zzj^) 

j¥=i 

The Macdonald-Koornwinder polynomials in N variables, P\{z), are 
labeled by partitions N parts, which parametrize the eigenvalue of 
Pq^t, he. the polynomials T’a(z) are characterized by the equations 

Vq^tPx{z) = dxPxiz) 

with 

N 

dx = J2 [q~^abcdt^'^-^-\q^^ - 1 ) + f-\q-^^ - 1 )] 

i=l 

together with the condition that the coefficient of in Pxiz) is 1. Of 
course the Macdonald-Koornwinder polynomials depend on the param¬ 
eter a, 6, c, d, t, q but we shall write this dependence explicitly (writing 
Px{z) = Pa(z; a, b, c, d; q; t)) only when needed. 

An important property of the Macdonald-Koornwinder polynomials 
(which could be used actually to give an alternative definition) is their 
orthogonality with respect to a certain scalar product. Let assume for 
convenience |g| < 1 and recall the usual notation for the q-Pochhammer 
symbol 

n-l j 

(104) (a;g)n = JJ(1 - ag*) (ai,a 2 ... ,aj;g)„ = JJ(aj;g)„. 

i=0 i=l 
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Let US define the kernel 


(105) Ap^{t,q-,z) 


TT 9)00 TT _ 9)00 _ 

ei,e2=±l 


and the scalar product on C[zf ,..., 


r ^ d 7 - 

{p, q) := JJ ^^A7v(t, g; z)p(z)g(z 

where the contour of integration C encircles the poles in ag^, bq^, cq^, dq^ 
{k G Z_|_) and excludes all the others. The Macdonald-Koornwinder 
polynomials corresponding to different partitions are orthogonal with 
respect to such scalar product 


(106) (Pa, = 0 if A 7 ^ p. 

When the partition A = {m}, i.e. when A consists of just a single part 
of length m, the Macdonald-Koornwinder polynomials are polynomials 
in one variable z independent of the variable t. They correspond (up to 
a factor and change of variables x = {z + z~^)/2) to the Askey-Wilson 
polynomials Pm(^; a, b, c, d|g)ll [ 2 ], [ 22 ] 

(107) -P{m}(z) OC Pm{z-, tt, b, c, djq). 


A.3. Askey-Wilson polynomials. Let us recall the dehnition of the 
Askey-Wilson polynomials in terms of a terminating hypergeometric 
function 


(108) 

Pn(z;a,b,c,djq) 


{ab, ac, ad; g). 


g ”, g” ^abcd, az, az 
ab, ac, ad 



where the basic hypergeometric functions r4>s are dehned by 
(109) 


bi,...,br 


;q,z = 


(Ul, . . . , Ur, g)n 

^ (g, 6 i,..., 6 ,;g), 


E 


;-i)”f (2 


■ 5+1—r 


As for the case of Macdonald-Koornwinder polynomials, whenever the 
parameters a, b, c, d are clear from the context we shall use a lighter 
notation writing Pn{z\q) = Pn{z; a, b, c, d\q). The Askey-Wilson polyno¬ 
mials satisfy orthogonality relations that correspond to the case N = 1 
of the orthogonality relation for Koornwinder-Macdonald polynomials. 


^In the literature the Askey-Wilson polynomials are usually considered as func¬ 


tion of a: = 


Z-\-Z 


2 
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Indeed for N = 1 the kernel Ai(f, g;z) becomes t independent and 
coincides with the Askey-Wilson kernel 


( 110 ) w{z]a,h,c,d\q) = 




{az, az~^, bz, bz~^, cz, cz~^, dz, dz~^] q)^ 


and one has 
( 111 ) 
where 
hn. 


dz 


Amz 


—w{z\q)pn{A(l)Pm{A(l) = hJn,m- 


(g" ^abcd] q)c 


(1 — g2>T.-ia5cd)(g”'+\ q^ab, q'^ac, q'^ad, q'^bc, q'^bd, q^cd] q)oo ' 


An important information we need in the paper is the asymptotic form 
of PmiA^) large values of the label m [221 Eq-7.5.13] 


( 112 ) pUAq) ~ z^B{z-A + z-^B{z) 


with 


B(z) 


{az, bz, cz, dz; g)oo 


The specialization c = d = 0 of the Askey-Wilson polynomials gives 
the so called Al-Salam Chihara polynomials [H [2S] 


(113) Qn{z; a, b\t) = Pn{z; a, b, 0, 0|t) 


The Al-Salam Chihara polynomials satisfy the following three terms 
recursion for 


(114) 

Qn+iiz) + ((n -|- b)t"' — 2z) Quiz) -|- (1 — t”')(l — abt^ ^)Qn-i{z) = 0. 


The recursion relations presented in Section |3]motivate the presentation 
of several contiguous/difference relations for Askey-Wilson polynomi¬ 
als, which are quite straightforward to prove. They are best stated 
using divided difference operators 


dx,yF{x, y) 


F{x,y) -F{y,x) 
{x-y) 


dxF{x) 


F{x) - F{x-A 
{x — X~A 
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and using the variable 2 ; instead of x 


2+Z~^ 

2 


Pm{z] a, b, c, d) = dc 4 [z] a, b, c, d)Pm{z] a, b, qc, d)) 

Pm-i{z; a, b, qc, qd) = dc,d a, b, c, d)Pm(z; a, b, qc, d)) 

Pm+i{z-,a,b,c,d) = dc 4 {u:^'^^^''^\z-,a,b,c,d)Pm{z-,a,b,qc,d)) 


Pm{z-, a, b, c, qd) = 
Pm+i{z-,a,b,c,d) = dz 
with 


^ 6 i(™’™)( 2 :; a, b, c, d)Pm{q^z; q^^a, q^b, q^^c, q^d)^ 
a, b, c, d)Pm{q^z; q^a, q^b, q^c, q^d)'^ 


uj^'^’^\z;a,b,c, d) 
a;("*-i’”^)(z;a, 6 ,c, d) 

u^^+^'^\z-,a,b,c, d) 
e^’^''^\z-,a,b,c, d) 
0 (’"+i’™)(z;a, 6 ,c, d) 


cd{z — c){zc — 1) 

(1 — cdq'^)zc 
1 

(1 - q'^){l - abq^-^) 

(1 — adq'^){l — bdq'^){z — c)(zc — 1) 
2; 

{abcq^ — z){l — az){l — bz){l — cz) 
gY(l — abq”^)(l — acq"^)(l — bcq'^)z‘^ 
(1 — az){l — bz){l — cz){l — dz) 

m ~ • 

q 2 2 :^ 


A.4. Cauchy formula for Macdonald-Koornwinder polynomi¬ 
als. We conclude this Appendix with the Cauchy formula for Macdonald- 
Koornwinder polynomials, due to Mimachi |30]. Let 

(115) n(z, x) = JJ {Zi + - Xj - xj^) 

l<i<N 

l<j<M 

and for a partition A C , i.e. a partition made of N parts such that 
Xi < M, define A, a partition of M parts given by 

Xj = 4h{AXi < N - i]. 

Theorem 7 ([30], Theorem 2.1). The Macdonald-Koornwinder poly¬ 
nomials P\{z) satisfy the equality 

(116) n(x,z)= ^ (-l)I^IPA(z;g;t)PA(x;t;g). 

\CM^ 

In the present paper we shall need just the case the case M = 1, 
which combined with the orthogonality property, allows to provide an 
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integral formula for the Macdonald-Koornwinder polynomials corre¬ 
sponding to a partition with a single column 


ri{N, m) 



N—m m 


involving the Askey-Wilson polynomials of basis t 
(117) 

ffi f (hoc 

= {{x + x)w{x\t)p^{x\t) 


where 

(118) 

{{x + X~^)'^ ,Prn) = 


{abcdq^'^] q) 


abq^, acq^, adq'^, bcq^, bdq^, cdq^] g)oo ’ 
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